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Abstract. The interesting but difficult phenomenology of supersymmetric models at the LHC and 
ILC demands a corresponding complexity and maturity from simulation tools. This includes multi- 
particle final states, reducible and irreducible backgrounds, spin correlations, real emission of pho- 
tons and gluons, virtual corrections etc. Most of these topics are included in the multi-particle 
Monte Carlo (MC) Event generators Madgraph, WHIZARD and Sherpa. A comparison of these 
codes is shown, with a special focus on the new release of WHIZARD. I show examples for the 
necessity of considering full matrix elements with all off-shell effects and interferences for multi- 
particle final states in supersymmetric models and give a status report on ongoing projects for 
simulations of SUSY processes at the LHC with these codes, including all of the abovementioned 
corrections. 

PACS. 11.30.Pb Supersymmetry - 12.38. Bx Perturbative Calculations - 12.60. Jv Supersymmetric 
models 



1 The need for multi-particle event 
generators 

Supersymmetry (SUSY) is the prime example for be- 
yond the Standard Model (SM) physics as a solution 
to the hierarchy problem. Compared to the SM, it 
doubles the particle spectrum, and the SUSY-breaking 
terms induce a vast number of new parameters. The 
analysis goal of the future experiments at the LHC and 
an ILC include mass measurements to determine the 
SUSY spectrum, the access of the new particles' spin 
encoded in angular correlations, and finally, by means 
of coupling measurements the proof that the new parti- 
cles really fill super-multiplets. Together with the com- 
plicated experimental set-up to dig out the informa- 
tions mentioned above, there is an utter need for pre- 
cise theoretical predictions of SUSY processes: first of 
all, to access them from the SM background, secondly 
since SUSY processes are themselves backgrounds for 
(more complicated) SUSY processes. Ultimately, the 
goal is to extract the SUSY Lagrangian parameters 
as precisely as possible to make predictions about the 
GUT structure and the SUSY breaking mechanism p] . 

In general, corrections to SUSY (or in general BSM) 
processes can be grouped into six categories: 1) loop 
corrections to SUSY production and decay processes, 
2) real gluon/photon radiation, 3) non-factorizable, 
maximally resonant photon exchange between produc- 
tion and decay, 4) off-shell kinematics for the signal 
process, 5) irreducible background from all other SUSY 
processes, and 6) reducible, but experimentally indis- 
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1. pr{b) and missing pr distribution for the process 
• fo&XiX? resulting from sbottom pair production at the 
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LHC. The SM background is in grey, the signal in blue. 



tinguishable background from SM processes. Topics 1) 
and 3) are dealt with in [2]. Here, we will focus on the 
last three points. 

SUSY processes reach a very high level of com- 



plexity already at tree level: e.g. the process e + e~ — > 
bbe + e~ XiXi ~ which is just X2X2 P a i r production - has 
~ 66, 500 Feynman diagrams. Several different produc- 
tion channels like XiX^ bib* and e^e* interfere with 
each other and need to be disentangled experimentally 
as well as in the simulation. One has to add SM back- 



grounds like e H 



And processes can be 
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much more complicated at the LHC, but also at the 
ILC. To make simulation of such complicated multi- 
particle processes feasible, among the following three 
levels of complexity usually one of the two first are 
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Fig. 2. Process: pp — » bbxiXu distributions of the harder 
and softer b jet in red and blue, respectively: on the left 
Pt, on the right in n. Full lines are full matrix elements, 
dashed ones Breit-Wigner approximation. 



used: first, there is the narrow-width approximation 
(NWA), where intermediate states are produced on- 
shell and multiplied by the corresponding branching 
ratio. Secondly, there is the Breit-Wigner approxima- 
tion (BWA), where the on-shell state's width is ac- 
counted for by folding in a Breit-Wigner propagator. 
Finally, there are the full matrix elements including all 
intermediate off-shell states contributing to the same 
exclusive final state. Programs like PYTHfA, HER- 
WIG, SUSYGEN and ISAJET use the first two ap- 
proaches. Here, we will show that this is not sufficient 
and that for BSM models full matrix elements have to 
used. 



Fig. 3. pt distributions of the four b jets in the process 
pp — + bbbbxiXu on the left: ordered according to their pr, 
on the right: ordered according to their centrality. 
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Fig. 4. px,b and t/>t distributions for the processes with 
two b jets in blue and with four b jets in red. 



2 SUSY simulations at the LHC 

For simulations of multi-particle SUSY processes at 
the LHC, the three packages WHIZARD 0HJ, Sherpa 
[5] and Madgraph [B] have been developed. All of these 
programs follow the SUSY Les Houches accord (SLHA) 
[TJ. Since the MSSM is a fairly complicated model with 
several thousand different vertices, it is mandatory to 
validate the correctness of its implementation in the 
codes. This has been done in [8] by means of unitarity 
and gauge invariance checks as well as a direct com- 
parison of the three programs. To test all phenomeno- 
logically relevant couplings, more than 500 different 
processes had to be checked; they are listed here [8j 
[5], which might serve as a standard reference. Fur- 
thermore, supersymmetric Ward- and Slavnov- Taylor 
identities have been checked [10]. 

To study the influence of interferences and off-shell 
effects as well as of radiative corrections, we chose a 
mSUGRA-inspired parameter point with non-universal 
scalar masses (note that the following discussion does 
by no means depend on that special point). The sbot- 
tom masses are 295 and 400 GeV, respectively, the 
lightest Higgs is directly above the LEP exclusion limit. 
It has a large (47 %) branching ratio of invisible decays 
to the neutralino LSP. The first generation squarks 
and sleptons are at 430 and 205 GeV, respectively. The 
point is compatible with all low-energy constraints like 



b — > S7, B s — > , Ap, g^ — 2, and cold dark mat- 

ter. The focus here lies on sbottom production with a 
BR(6i -> bxl) = 43.2 %. 

Fig.[T]shows the parton-level distributions for pp — > 
bbxiXi> where we used standard cuts of pr,b > 20 GeV, 
r\h < 4, and AR^ _> 0.4. The main SM background 
comes from pp — > bbvv and is shown in gray. The sig- 
nal is in blue, the pT,b distribution on the left, the 
j>T distribution on the right. Here, the bumb coming 
from the sbottom production is clearly visible. Gener- 
ally, the signal jets are harder than the jets from the 
SM background. In Fig. [2J the px.b and % distribu- 
tions are plotted for the harder (red) and the softer 
(blue) of the two jets. From the right plot we see that 
the harder jet is much more central as expected from 
phase space restrictions. The full lines here denote us- 
age of full matrix elements, while the dashed lines are 
the BWA. There is a clear discrepancy between the 
curves, stemming mainly from bbZ* diagrams; these 
discrepancies are only in the low-px region which is 
anyhow cut out. The reason why off-shell decays are 
irrelevant here is the large mass splitting from 295 to 
45 GeV between the sbottom and the LSP. For gen- 
eral SUSY decay cascades this is not the case as will 
be explained in section 01 

As a next step, we studied the influence of real cor- 
rections on the detection of the signal, he. bottom-jet 
radiation from the initial state (g — > bb splitting) as 
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Fig. 5. Top: energy of the b jet in the process e + e~ — > 
bb + at the ILC. In red the NWA, in light gray the SM 
background, and in dark gray the signal with full matrix 
elements on top of it. ISR and beamstrahlung are included. 
Bottom: same, but f£ distribution. 

combinatorial background. Using the above cuts for 
the b jets, the perturbation series is well-behaved, the 
cross section ceases from 1177 fb to 130.7 fb when in- 
cluding the two bottom jets. The process with four 
b jets, pp — > bbbbxiXii is a t the border of feasibility 
for SUSY LHC simulations with 32,000 diagrams, 22 
color flows and several thousand phase space channels. 
The simulation has been performed with WHIZARD, 
which is well-suited for physics beyond the SM pTj . 
Fig. [3] shows the pr distributions of the four b jets, 
ordered according to their px on the left and to their 
pseudorapidity on the right. The right figure shows 
that only the most forward jet is considerably softer, 
and that a forward discrimination between ISR and 
signal (decay) jets can be quite intricate. Fig. [3J on 
the other hand, shows that the jet structure (left hand 
side) itself does not change considerably from the case 
with two (blue) and four b jets (red). Since more hard 
jets have to be produced, the PDFs decrease the max- 
imum a little bit. The missing px is shifted to lower 
values because light particles balance out the event. 



3 SUSY simulations at the ILC 

At the ILC, we consider the same final state, e + e~ — > 
bbxiXi- Here, the production is electroweak, whence 



a lot more of intermediate states contribute: Zh, ZH, 

HA, x?x°, Xixi XiXl hk, bJ$- Especially, the back- 
ground from heavy Higgses and the heavier neutrali- 
nos is quite severe. The irreducible SM background is 
mainly from WW fusion: e + e~ — > bbvv. To extract 
the sbottom from the SUSY backgrounds, one has to 
cut out the regions in the M b i spectrum from 150 to 
250 GeV as well as from 350 to 800 GeV. Applying 
these cuts, the result for the cross section of the signal 
is 0.487 fb (0.375 fb with ISR and beamstrahlung); the 
NWA gives 2.314 fb, which is one order of magnitude 
away from the full result. The reason is that interfer- 
ences with heavy neutralino three-body decays affect 
the decay kinematics of the sbottoms. Hence, for ILC 
it is absolutely mandatory to use full matrix elements 
and multi-particle event generators for SUSY simula- 
tions. 



4 Origin of off-shell effects 

In this section, we briefly explain the reason for the 
appearance of large off-shell and interference effects 
for SUSY processes at the LHC (for a clear deriva- 
tion of these effects, see [32]). Although, here we fo- 
cus on SUSY, all new physics scenarios which have 
a conserved matter parity (responsible for dark mat- 
ter) and a complicated spectrum without too large hi- 
erarchies in common, share the same features. Their 
phenomenology grossly consists of decay chains with 
quasi-degeneracies among mother and daughter parti- 
cles. Including the full momentum dependence of the 
decay matrix elements affects the momentum depen- 
dence of the propagators: performing phase space in- 
tegrals does not result in simple Breit-Wigner terms 
that factor out. Although naively, the importance of 
off-shell effects is expected to scale with the width to 
mass ratio FjM for the resonances, these factors can 
be largely enhanced by powers of the (small) inverse 
velocity of the daughters, l//3 n , coming from the near- 
degeneracies of mothers and daughters. Furthermore, 
interferences with diagrams where only one of the two 
decay chains is resonant can be quite large. 

Altogether, these effects lead to drastic deviations 
from the NWA: effective branching ratios (which are 
measured via their exclusive final states) are shifted 
by 20 — 40 %. Charge and chirality asymmetries could 
appear in the decays of gluinos: e.g. gg — > bb\bb\ vs. 
bbibbl could have different effective branching ratios 
since they have different subsequent decay chains which 
have different interferences with non-resonant diagrams 
A chirality asymmetry is a non-equality between the 
two decays gg qL<lLjj vs. qnqujj for the same rea- 
sons as above. 

In Fig. [S] (taken from [T2]) the effect of using full 
matrix elements are shown for the process ud — ► Xi9 
with subsequent decays of the gluino. The left plot 
shows the deviations of the effective branching ratios 
for the gluino measured with the help of exclusive fi- 
nal states as a function of the ratio of the squark over 
the gluino mass. As this ratio approaches one, off-shell 
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Fig. 6. Off-shell and interference effects for SUSY decay 
chains at the LHC. Left: shifts to the gluino BR as a func- 
tion of rriq/mg. Right: Decay asymmetry (g — > q~Lq)/(g — > 
q~B.q); from [12] • The shaded blue area covers all SPS points. 



and interference effects give corrections of the order of 
20 — 40 % compared to the NWA. This near-degenerate 
squark-gluino parameter region is quite generic, since 
all SPS-standard candle points [13] lie in the blue band. 
So for all of these points, there are large deviations 
from the NWA. The right plot shows a chirality asym- 
metry, i.e. an asymmetry between the decays g — > qj,q 
and g -> q R q. 

In [14] , all of the off-shell and interference effects 
discussed briefly above, are studied systematically for 
all kinds of different SUSY decay chains and observ- 
ables. Also, the size of these effects as a function of the 
region of parameter space is investigated. 



5 Conclusions 

In conclusion, we validated the three multi-particle 
event generators for MSSM processes, Madgraph, Shcr- 
pa, and WHIZARD, to produce correct results. We 
studied the influence of using full matrix elements for 
exclusive four- and six-particle final states in SUSY 
production processes at the LHC and ILC compared 
to the narrow-width or Breit-Wigner approximation. If 
there are several different diagram classes that can be- 
come singly- or doubly-resonant, interferences of these 
classes result in deviations from the NWA up to an 
order of magnitude. Furthermore, off-shell effects are 
important for nearly-degenerate mother-daughter con- 
stellations, as is typical in SUSY decay chains (but also 
in all BSM scenarios with a conserved parity and sim- 
ilar spectra). Especially, when cuts are mandatory to 
disentangle signals and backgrounds, interference and 
off-shell effects have to be studied by multi-particle 
event generators. 
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